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ABSTRACT 

The Ultra Luminous InfraRed Galaxy (ULIRG) Mrk231 reveals up to seven rotational lines of water (H2O) in emission, including a very high- 
lying (iiuppcr = 640 K) line detected at a Aa level, within the Herschel/SPIRE wavelength range (190 < /l(//m) < 640), whereas PACS observations 
show one H2O line at 78 yum in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused 
by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, 
which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. 
Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate 
the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 
300 - 640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius= 1 10 - 180 
pc) and warm {T^^st = 85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, ^(1120) ~ 5 x 10" cm"^, 
found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" 
where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, 
could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 jira line observed with PACS 
shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer 
et al., this volume). Additional PACS/HIFI observations of H2O lines are required to constrain the kinematics of the nuclear component, as well 
as the distribution of H2O relative to the warm dust. 

Key words. Line: formation - Galaxies: ISM - Galaxies: individual: Markarian 231 ~ Infrared: galaxies - Submillimeter: galaxies 



1. Introduction 

One key question in the study of composite infrared (IR) merg- 
ing galaxies and quasi-stellar objects (QSOs) is what fraction of 
their luminosity is generated in the nuclear region (< 200 pc) 
associated with the Active Galactic Nucleus (AGN) and a pos- 
sible extreme nuclear starburst, and what fraction arises from 
a more extended kp c-scale starburst (e.g. lArmus et al.L l2007t 
IVeilleuxetalT.l200S<) . The ULIRG Markarian 231 (Mrk 231) is 
the most luminous (L ~ 4x 10'^ L©) galaxy in the local Universe 
iz < 0.1), and thus provides a unique template for such stud- 
ies. Since the bulk of the luminosity in ULIRGs arises at far-IR 
wavelengths, where sub-arc-second resolution observations are 
not available, an alternative technique is required to constrain the 
compactness of the far-IR emission and its physical origin. 

In a previous work bas ed on observations with the Inf rared 
Space Observatory (ISO), iGonzalez- Alfonso et al.l (l2008l here- 
after G-A08) have argued that the observation of molecular 
species such as OH and H2O at far-IR wavelengths is ideal for 
such a purpose, because their high-lying rotational levels are 
pumped through absorption of far-IR radiation and the observ- 
able excitation is then sensitive to the far-IR radiation density 
that in turn depends on the compactness of the far-IR contin- 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



uum source. In addition, these molecular observations shed light 
on the dominant chemistry in those nuclear regions. G-A08 re- 
ported the ISO detection of 3 high-lying H2O lines, relevant 
upper limits over the entire ISO spectrum, and also high-lying 
OH lines, indicating the ocurrence of a compact-luminous far- 
IR component. 

With their high sensitivity, spectral resolution, and wave- 
length coverage, the Herschel (Pilbratt et al., 2010) instruments 
are ideal for extending our previous study to additional key 
lines in the far-IR /submillimeter. As part of the HerCULES Key 
Programme (see Ivan der Werf et all 1201 Ol this volume, here- 
afte rvdWlO) , we report in this Letter the Herschel SPIRE/PACS 
(Griffin et all 1201 Ot iPoditsch etalL [2010) detection and first 
analysis of several H2O lines in Mrk 231, which supports the 
conclusions of G-A08 and gives additional clues to the origin of 
H2O in this ULIRG. We adopt a distance to Mrk 23 1 of 192 Mpc 
{Ho = 70 km s-i Mpc"', Qa = 0.73, and z = 0.04217). 



2. Observations 

The SPIRE FTS observations of Mrk 231 (vdWlO) were con- 
ducted on December 9th, 2009. The PACS observation of one 
H2O line was conducted on November 8th, 2009, as part of the 
SHINING Key Programme, and kindly provided for the present 
study. Details on data re duction, calibra t ion, an d line extraction 
are given in vdWlO and lSwinvard et all (l2010h . Excerpts of the 
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Fig. 1. Comparison between the observed spectra (black/green 
histograms: unapodized/apodized spectrum, see vdWlO) and re- 
sults for the reference model discussed in §3 (dark blue lines; a 
Gaussian instrumental line shape with FWHM - 0.048 cm ' is 
used for simplicity). The 312221 line, shown in the same panel 
as the 3"21312 line, is blended with '^CO (10-9) (vdWlO). The 
red segment in each panel indicates the FWHM of an unresolved 
line. The velocity scale has been calculated with respect to the 
systemic redshift of z = 0.04217. 



spectrum around the H2O lines are displayed in Fig.[Tl and an en- 
ergy level diagram indicating the lines detected with SPIRE, the 
one detected with PACS, and those detected with ISO is shown 
in Fig.|2] Table [1] lists the line fluxes. Figure [1] also shows the 
results of our reference model, discussed below. 

H^^^O could be marginaUy detected at 250.0 yum (220211) 
(see vdWlO); however, this identification should be confirmed as 
the feature is shifted by 150 km s"' from the nominal line wave- 
length. Confirmation of Hl^^O would be important as its pres- 
ence would support the strong enhan cement of "^O i n Mrk 23 1 
derived from "*OH observations (Fis cher et al.Ll2010l this issue, 
hereafter FIO). A number of the H2O lines in Fig. [T]are blended 
with the C'^^O fines 9-8 (303.57 jjm), 10-9 (273.24 jum), 
and 11-10 (248.43 //m). Contamination by C'^O is minimal, 
however, as the lower-lying 6 - 5, 7 - 6, and 8-7 lines are not 
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Fig. 2. Energy level diagram for H2O, showing the de- 
tected/undetected (blue arrows/lines) lines with SPIRE, the line 
detected with PACS (green) and those detected by ISO (light 
blue). Dashed red arrows indicate the main pumping paths for 
the high-lying lines observed with SPIRE. Upward (downward) 
arrows: absorption (emission) lines. 



detected. PACS observations show a broad absorption feature at 
121.7 ;um, nearly coincident with the H2O 432423 line (FIO). 
However, this feature is probably contaminated by HE (2 - 1) at 
the same wavelength, as the 1-0 line is detected with SPIRE 
(vdWlO); HE has been prev iously detected towar ds the Galactic 
star forming region Sgr B2 (Neu feld et aU[l997l) . Therefore the 
121.7 jum feature is not used for the H2O analysis below. 

Of particular interest is the detection of the very high-lying 
5235 14 H2O line at a Act level, which we have verified by repro- 
cessing the data with no correction applied for the instrument 
spectral efficiency and by comparing these data with reduced 
observations of the dark sky over the same spectral range. The 
ground state lines (P-H2O 111000 in Fig. [Hand 0-H2O 110101) 
are not detected. The lines detected with SPIRE are all in emis- 
sion and peak at central velocities, in contrast to the low-lying 
OH lines that show P-Cygni profiles characteristic of an extreme 
molecular outflow (FIO). The red horizontal lines indicate the 
FWHM of an unresolved line, and show that the H2O lines de- 
tected with SPIRE are barely resolved. The 0-H2O 423312 line 
detected with the higher spectral resolution of PACS is in ab- 
sorption and well resolved, showing a central main body and, 
apparently, a relatively weak blue-shifted wing extending up to 
-800 km s ' and possible low level emission from the receding 
gas. These wings should be confirmed with additional observa- 
tions of H2O absorption lines, as the limited wavelength cov- 
erage of the 423312 line makes the adopted baseline uncertain. 
The detection would indicate that H2O also participates in the 
prominent outflow detected in OH (FIO). 

While the shape of the 4233 12 absorption line shows the cen- 
troid of the main body slightly blue-shifted (by ~ -70 km s"'), 
some lines observed in emission tend to show, on the contrary, 
a slight red-shift of their centroid (up to 100 km s '). This ef- 
fect could be related to systematic motions, i.e. a low-velocity 
nuclear-scale outflow, and will be explored in the future with ad- 
ditional high spectral resolution observations. 



3. Analysis 

The observed pattern of line emission cannot be explained in 
terms of pure colfisional e xcitation (G-A08). Adopting the colli- 
sional rates of iFaure et aP (I2007.) with Tk > 200 1<3 and «(H2) = 



An LTE ortho-to-para H2 ratio according to Tj is assumed. 
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Table 1. Observed and modeled line fluxes. 



Line 



A Flux" Model 

(yum) (Jy km s"') (Jy km s"') 



P-H2O 111000 


269.27 


<300 


177 


0-H2O 110101 


538.29 


<400 


146 


P-H2O 202111 


303.46 


718(110) 


660 


P-H2O 211202 


398.64 


415 (43) 


365 


P-H2O 220211 


243.97 


342 (92) 


372 


0-H2O 312303 


273.19 


400 (130) 


438 


0-H2O 321312 


257.79 


631 (47) 


657 


P-H2O 422413 


248.25 


361 (38) 


381 


0-H2O 523514 


212.53 


287 (83) 


269 


0-H2O 423312 


78.74 


-910'' (60) 


-664 



" Numbers in parenthesis indicate the estimated uncertainties. 
'' Includes the absorption in the high-velocity blue-shifted wing, 
which accounts for 220 Jy km s"'. 



10 - 10 cm , and ignoring radiative pumping, the models 
that account within a factor of 2 for the high-lying 523514 and 
422413 lines predict fluxes for the low-lying lines that exceed 
the observed values by factors of > 10. Thus the observed line 
ratios indicate an excitation mechanism that favours the emission 
in the high-lying (> 300 K) lines at the expense of the low-lying 
lines. This does not imply that the H2O lines are not formed in 
warm/dense regions, but just that the dominant excitation mech- 
anism for the high-lying lines is not collisional. 

Such a mechanism is the pumping through absorption 
of dust-emitted far-IR photons, which efficiently pumps the 
high-lying 3"21312/422413/5'23514Hnes through absorptions at 
75.4/57.6/45.1 fim in the strong 3"21212/422313/5'23414 lines 
(Fig. |2]), of which the lower backbone levels are preferently pop- 
ulated. This requires a strong continuum component at 30 - 70 
fim. However, this component cannot dominate the emission at 
> 130 |/m, as it would produce strong H2O absorptions in that 
wavelength range that are not observed (G-A08). The data then 
support the occurrence of both a warm/compact component with 
moderate opacity, and a second colder component naturally as- 
sociated with the more extended 1-kpc starburst that dominates 
the emission at long wavelengths. Our proposed SED decompo- 
sition is shown in Fig. [3^, and defines the reference model {M^^i) 
with parameters as listed in Table |2| (0 a hot component {He) 
with Tdust = 150 - 400 K dominates the emission at /I < 20 
//m; (//) a warm (95 K) and compact (R - 120 pc) component 
iWc) dominates at 20jum < A < lOfim; (Hi) an extended 1-kpc 
component (Ec), with Tjust = 40 K, accounts for most of the 
continuum at /I > 70 fj.m. The Wc, with Lw ~ 1.9 x 10'^ Lq, is 
responsible for the observed high-lying H2O line emission. 

Calculati ons for H7O were carried out using the code 
described in iGonzalez- Alfonso & Cernicharol ( Il999h . In M^^f 
(Table |2]i, line broadening is caused by microturbulence. We 
have adopted a "mixed" approach (i.e. the H2O molecules are 
evenly mixed with dust, G-A08), discussed below. An ortho-to- 
para H2O abundance ratio of 3 is assumed. For those H2O lines 
observed in emission, Fig.[3]3 compares the expected fluxes from 
the Wc (in violet) with the observed fluxes (in red). Collisional 
excitation is included with gas at T^ = 150 K and n(H2) - 
1 .5 X lO*" cm"^ but, even for these shock-like condition^ it has a 
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Fig. 3. Comparison between observations and results for the ref- 
erence m odel, a) Continuum emission from Mrk 231. Spitzer 
IRS data dArmus et aI.L |2007|) . ISO data (G-A08), and SPIRE 
data (red spectrum) are shown. Flux densitie s at 800 and 1100 
jum are taken from Ro che & Chandlei (1 19931 corrected for non- 
thermal emission), at 450 yum from'Rigopoulou et al.l (11996 '). and 
at 350 yum from Yang & Phillips (2007). b) H2O emission. The 
high-lying lines are generated in the warm/compact region (vio- 
let), whereas the extended kpc-scale starburst (green) is expected 
to contribute to the low-lying lines. 



low efifect on the H2O level populations and line fluxes as these 
are mostly determined by the strong radiation field. The high- 
lying lines are reproduced with the Wc, but there is a model 
deficit of emission from the low-lying lines. This deficit indi- 
cates that the Ec contributes to those low-lying lines (green); 
both radiative and collisional excitation, the latter significant for 
£^upper ^ 200 K, are included in the model. 

The high-lying 422413 and 523514 lines are (nearly) opti- 
cally thin in the Wc, so that their expected fluxes are sensitive 
to A^(H20). Table|2]shows that, despite the strong radiation field 
in this region, a high A^(H20) is required to account for the ob- 
served fluxes. Assuming a gas-to-dust mass ratio of 100, the av- 
erage H2O abundance relative to H2 in this Wc is 0.8 x 10"^. 

The H2O line at 78.7 /im observed by PACS is consis- 
tently predicted in absorption, and its flux is reasonably repro- 
duced (Table[T]i, given that the blue-shifted wing is not modeled. 
However, the observed line shape suggests systematic motions 
that are not included in our model. In general, we expect that 
lines with A < 120 //m are observed in absorption, while lines 
with A > 120 jjm are in emission (with exceptions due to the 
pumping details and level energies). 

Concerning the ISO lines (G-A08), M,ef also matches the 
absorption in the 2201 1 1 line at 101 /im, but underestimates the 
absorption in the 330221 and 331220 lines at 66.4 and 67.1 yum 
by a factor of 2. As mentioned above, our model uses a mixed 
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Table 2. Parameters of the reference model (Mref). 



Components 


Hot" (He) 


Warm (Wc) 


Extended (Ec) 


Radius (pc) 


23 


120 


610 


Tdmt (K) 


400 - 150 


95 


41 


TlOO/mi 


0.4 


1.0 


0.5 


i(Lo) 


7.5 X 10" 


1.9x10'^ 


9.6x10" 


Gas Mass*" (Mq) 


1.9 X 10^ 


5.9 X 10** 


7.7 X 10' 


A'(H20) (cm--) 


- 


5.2x10" 


2.0 X 10'" 


V,„A (km s-') 


- 


60" 


40 


n(ll2f (cm-3) 


- 


1.5x10" 


5 X 10"^ 


71. (K) 


- 


150 


100 



" The hot component does not contribute to the H2O emission. 

*■ A gas-to-dust mass ratio of 100 is assumed. 

■^ These parameters are not well determined for Wc, as the excitation 

is dominated by radiative pumping. 
" From.Do wnes & S olomon, Cl998h . 



approach, which implies that for given values of A^(H20), the 
absorbing H2O lines are relatively weak, as molecules located 
deep inside the source do not contribute to the absorption fea- 
tures. Conversely, if a screen approach is adopted (i.e. a H2O 
shell surrounds the continuum source), the absorption lines be- 
come much stronger In both approaches, the A^(H20) required 
to match the lines observed with SPIRE are similar Thus an 
analysis that combines absorption and emission lines is a pow- 
erful tool to establish the distribution of H2O relative to the 
warm dust responsible for the excitation. The screen version of 
Mref yields absorption in 330221, 331220, and other transitions 
that overestimate the observed values or upper limits. Therefore, 
our preliminary result is that a combination of both the mixed 
and screen scenarios best describes that observed data, with the 
mixed version favoured. Nevertheless, it remains unclear what 
fraction of the observed 330221 and 331220 absorption arises 
from the outflow detected in OH (FIO). 

By increasing Re, the radius of Wc, to 170 pc, and keeping 
Lw constant, the 321312/4233121ine strengths become overesti- 
mated by 50/25%, and the 211202/5'23514 intensities are un- 
derestimated by 30%; given the simplicity of our spherically 
symmetric models, we estimate a size for Wc in the range 
Re = 1 10 - 180 pc. A lower limit for the luminosity arising from 
Wc is estimated by decreasing Tdust to 85 - 80 K and increas- 
ing A^(H20) to > 10'** cm"^, which results in too weak 423312 
and ISO absorption lines. We estimate that the mid- and far-IR 
emissions from the nuclear region account for more than 45% of 
the bolometric luminosity; observations at 60 - 200 fim are re- 
quired to better constrain that value and to establish a firm upper 
limit. Results are more uncertain for Ec, and we may expect that 
its contribution to the H2O emission arises from its innermost 
region. 

4. Discussion 

The extreme nature of the nuclear region in Mrk 23 1 is well illus- 
trate d by compari n g its S PIRE spectrum with that of the Orion 
Bar dHabart et al.L l2010l this issue), the prototypical Galactic 
PDR. The Orion Bar spectrum shows CO lines a factor of > 50 
stronger than the H2O lines, while in Mrk 23 1 the H2O and CO 
lines have comparable strengths (vdWlO). This contrast will be 
still higher in the nuclear region, provided that a significant frac- 
tion of the CO emission in Mrk 23 1 arises from a more extended 
region. Thus the H20-to-CO line intensity ratios in the SPIRE 



wavelength range are an excellent diagnostic of extragalactic 
compact/warm far-IR continuum sources with unusually high 
amounts of H2O. 

The above comparison also indicates that the nuclear region 
of Mrk 23 1 cannot be interpreted as an ensemble of classical 
PDRs. Three main scenarios are proposed to explain such high 
amounts of H2O: (/) widespread shocks/cosmic rays: although 
the H2O lines peak around the systemic velocity, outflows of 
~ 100 km s' are not ruled out by our data, and indeed some 
indications in the H2O line shapes of systematic motions have 
been found; an enhanced cosmic ray flux could also have an im- 
portant impact on the nuclear chemistry. (//) XDR chemistry: 
our derived H2O abundanc e of ~ 10'^ is in very good agreement 
with XDR model results bv lMeiierink & Sp aans ( 2005, their Fig. 
3, Model 3), as well as with our preliminary estimate of the H2O 
spatial distribution. (///) An undepleted chemistry, where H2O 
that formed on grain mantles is released into the gas phase, as in 
Galactic hot cores; in support of this scenario, the derived Tdust 
in Wc is close to the evaporation temperature of solid H2O. All 
three scenarios are probably taking place, and the identification 
of the dominant process requires a multi-species analysis. 

The nuclear region traced by the high-lying H2O lines has 
a size similar to the nuclear disk (or outflow' ) observ ed at radio 
wavelengths and H I 21 cm by Carilli et al.l (Il998l their Figs. 
3 & 7), suggesting a close physical correspondence. From He 
and Wc, the nuclear surface brightness (~ 1.5 X lO'^ L© kpc"^) 
exceeds th e highest value s attain ed in starburst on sp atial scales 
> 100 pc (M eureretall 119971: iDavies et aU 120071) . while the 
nuclear luminosity-to-mass ratio (L/M ~ 4 x 10^ Lq/M q) ex- 
ceeds the limi t for a starburst estim ated bv lScovilIa(l2003h . From 
near-IR data. iDavies et al.l (l2007h estimated a starburst luminos- 
ity from a similarly sized region of < 7 x lO" L©; according to 
the joint luminosity of our He and Wc, the AGN would account 
for at least 50% of the output power in Mrk 23 1 . 
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